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ABSTRACT 

The total amplification of a source inside a caustic curve of a binary lens is 
no less than 3. Here we show that the infimum amplification 3 is satisfied by a 
family of binary lenses where the source position is at the mid-point between the 
lens positions independently of the mass ratio which parameterizes the family. 
We present a new proof of an underlying constraint that the total amplification 
of the two positive images is bigger than that of the three negative images by one 
inside a caustic. We show that a similar constraint holds for an arbitrary class of 
n-point lens systems for the sources in the 'maximal domains'. We introduce the 
notion that a source plane consists of graded caustic domains and the 'maximal 
domain' is the area of the source plane where a source star results in the maximum 
n 2 + 1 images. We show that the infimum amplification of a three point lens is 
7, and it is bigger than n 2 + 1 — n for n > 4. 

subject headings: gravitational lensing - planetary systems - stars: binaries: 
general 
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1. Introduction 



Gravitational lensing has been known for three quarters of a century and has 
been extensively used in studies of quasars and clusters of galaxies. In these cases, 
the lenses are compound lenses consisting of extended objects such as galaxies and 
clusters of galaxies. The recent gravitational microlensing experiments looking 
for baryonic dark matter (Alcock, et al., 1993, Aubourg, et al., 1993, Udalski, et 
al., 1993) study nearby stars 50kpc) as source stars, however, and the lenses 
are point lenses. They are mostly single lenses but some of them have turned 
out to be binary lenses (Udalski, et al, 1994, Bennett, et al, 1995, Alcock, et 
al, 1995, Alard, et al, 1995). The detection rate of binary lenses is expected 
to be about 10% according to Mao and Paczyhski (1991). One interesting class 
of binary lenses is that of planetary systems, and there is an effort to look for 
extra solar earth mass planets using gravitational microlensing (Tytler et al, 
1995). Obviously, the observational advantage of having to deal with point lens 
systems is that the simplicity allows high precision experiments. From a theorist's 
point of view, the current and proposed microlensing experiments have brought 
a necessity to study fine details of the cleanest lens systems, which in turn would 
lead to a more sophisticated understanding of more uncertain lens systems. 

It is well known that when two masses are nearby, they produce a gravita- 
tional lens with a rich structure that depends on the mass ratio and separation 
between them and the lens structure is best represented by the caustics. The 
caustic curves are closed cuspy loops, and neither do they self-intersect nor nest 
each other. (If n > 3, the lenses are general enough to feature self-intersections 
and nesting of the caustic loops. See figure 1 for an example of a ternary lens.) A 
caustic curve of a binary lens defines an inside and an outside and the number of 
images of a source is five inside a caustic loop and three outside. Infinity lies in the 
outside domain and the number of the images of a source at infinity is three, one 
at the source position and the other two at the lens points. Obviously, the source 
at infinity is unamplified and the total microlensing amplification Atot(oo) = 1. 
The amplification of the image at the source position is one and that of the im- 
age at each lens position is zero. The amplification of a positive image (= an 
image with positive parity) is always no less than 1, and the total amplification 
Atot > 1 in the outside domain. In an inside domain, a source produces two 
positive images and hence A to t > 2. In fact, A to t > 3 inside a caustic as shown 
recently by Witt and Mao (1995). Their main result was the following relation. 




(1.1) 
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where J is the Jacobian determinant of the lens equation. In order to demonstrate 
(1.1), the authors converted the binary lens equation inside a caustic into a 
fifth order polynomial equation of J -1 using an algebraic computational package 
Mathematica. They also made an empirical suggestion that there is only one 
binary lens configuration that accommodates the infimum amplification 3 inside 
a caustic. 

Here we show that there are actually one-parameter family of binary lenses 
where A tot > 3 is saturated inside a caustic. We also present a simple algebraic 
proof of (1.1) and show that a similar "sum rule" holds in the "maximal domains" 
of arbitrary n-point lens systems. This leads to a necessity to look at the source 
planes as consisting of graded caustic domains {V m }. In section 4, we discuss 
the hierarchical structure of the source plane of an arbitrary n-point lens and 
infimum microlensing amplifications of the sources in the "maximal domains." 



2. Inside the Caustics of Binary Lenses 

A lens equation is a mapping from an image position to a source position 
on the lens plane transverse to the line of sight. If we set the distance scale of 
the lens plane by the Einstein ring radius of the total mass, the normalized lens 
equation is given by (Bourassa, Kantowski and Norton, 1973) 

U = Z-— --— r , (2.1) 

(Z-Xl) (Z-X2) 

where z and u are the image and source positions in complex coordinates, x\ 
and X2 are the lens positions on the real axis, and e\ and €2 (ei + 62 = 1) are 
the fractional masses of the lens elements. The Jacobian matrix elements of the 
lens equation (2.1) are 

d z u = 1 , dzu = 61 » + 62 » = R , (2.2) 

\Z — X\Y \Z — X2) 

and the Jacobian determinant is J = |dzcu| 2 — |<9zu;| 2 = 1 — |<9^| 2 < 1. Lensing is 
an inverse mapping of the lens equation and the inverse Jacobian matrix elements 
are given by 

1 K 

duZ = — , d z = -- . (2.3) 

For a positive image, < J < 1 and the amplification A=|J| _1 >1. A source 
inside a caustic has two positive images, and the total amplification of the positive 
images A + > 2. Therefore, combined with (1.1), the total amplification A to t > 3. 
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In order to find all the lens parameters where the inequality At Q t > 3 is 
saturated, we note that the minimum Atot = 3 can occur only if both of the 
positive images have J = 1 (or d^oJ = 0). They are the two finite limit points of 
a binary lens (Rhie, 1995a). If I = \x\ — X2\ denotes the separation between the 
lens masses, 

Z*± = £acm ± V e l e 2 I 5 X acm = 6iX 2 + € 2 Xi . (2.4) 

The limit points z*± are the images of u*± which we can calculate from the lens 
equation (2.1). If we let e = e\ — e 2 , 

w*± = £acm - ^(x 2 - xi) ± iy/eiei (l ~ j^j , (2-5) 

and u;* + and cu*_ have to be degenerate in order to be one source producing two 
images at the limit points. 



I = V% ; w* = -(#1 + x 2 ) . (2.6) 



In other words, a family of binary lenses with seperation / = \/2 accommodate a 
source position with the infimum amplification 3 inside a caustic and the source 
position is at the mid-point of the lens positions irrespectively of the mass ratio 
that parameterizes the family. 

The fact that there are two positive images implies that the source in (2.6) 
is inside a caustic. (In a binary lensing, the number of negative images is always 
bigger than that of positive images by one.) The three negative images are on 
the lens axis and all have J = —3. If the source is at the origin [x\ + x 2 = 0), 
the lens equation for the negative images (on the real axis) is given by 

= x(x 2 - \) + ^= , (2.7) 

where x 2 > x\ is assumed without loss of generality. The Jacobian determinant 
is given by 

J = 1 - ( ( ?~ a T2 2 + Cl ^n 2 - 1 " N 2 • (2-8) 



(x — xi) 2 (x — x 2 ) 2 

When the lens is symmetric (e = 0), the images are at x = , ± a/3/2 and it is 
easy to verify that J = — 3 for each image. For an arbitrary e, one can reduce the 
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numerator and denominator of \k\ 2 into second order polynomials using (2.7) 
to find that \k\ 2 = 4 nd so J = —3. The family of lenses defined by / = \[2 
have caustics consisting of either one connected loop or two. (For a three-loop 
caustic, / < 1. See Rhie, 1995a). Two-loop caustics occur when the minor mass is 
small, and the source position for the infimum 3 is in the caustic near the minor 
mass. Inside the caustic (very) near the dominant mass, the amplification goes 
to infinity as the mass asymmetry grows because the caustic approaches that of 
a single lens. 

Here we demonstrate that (1.1) does not hold when the source is outside - 
the non-maximal domain of a binary lens - except when the source is at infinity. 
Let's consider a source located at the mid-point of a symmetric lens. When / < 2, 
the source lies inside the caustic, and there are five images: two off the axis 
(z = ±Vd 2 — 1; d = 1/2), and the other three (z = 0, ±V d 2 + 1) on the axis. As 
the separation / grows beyond the bifurcation condition 1 = 2, the caustic splits 
into two closed loops and the source falls in the outside domain. Thus, two out of 
the five images should disappear, and they are the two solutions z = ±V d 2 — 1, 
which become real at / = 2. Both have J = — 4d 2 (d 2 — 1) (< OforZ > 2). Therefore, 
the other three images have ^images J 1 > 1 ( — * 1 as <i — > oo). Of course, the 
total amplification of the three images of a source outside can be dominated by 
negative images and ^images J 1 < 1- An example is a source near the off-axis 
cusps which appear when the caustic bifurcates into three loops. 



3. ^2 J 1 = 1 for the Maximum Number of Images 

The lens equation (2.1) is a mapping from an image position to a source 
position: z,z i— > uj, and lensing is the inverse of the lens mapping: ui, u i— > z 
which is either triple- valued or quintuple- valued. If we let z\ = z — x\, zi = z — 
and z = z — x &cm , the binary lens equation looks quite simple. 

u = z- — = z-f(z). (3.1) 

We note that the effect of the masses ej : j = 1, 2 is coupled only to the conjugate 
variable z and the relation between uj and z is universal. We will see below that 
this universality is the underlying reason for the constraint (1.1), and thus, it 
holds not only for binary lenses but also for arbitrary n-point lenses. 

In order to find the inverse function z(u>, u>), we replace z in (3.1) by ui + f(z): 
uj = z — f{u) + f(z)) to find that z satisfies a fifth order polynomial equation, 
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= g(z; uj, uj) (Witt, 1990). If uj; j = o, 1, 2 is defined similarly to Zj, 

g(z; uj, uj) = (z - uj)(z + Q\Z\Z2){z + O2Z1Z2) - z\Z2 (z a + G) z\z<i) . (3.2) 

g = is an analytic equation in z and always has five solutions. Thus, g = 
is not equivalent to the original lens equation when the source is outside where 
there are only three images. However, in the inside domains we are interested 
in here, the sources result in the maximum number of images 5 and so we can 
consider = g as the lens equation. If aj(uj, uj); j = 1, 5 are the zeros of g, we can 
factorize g as g(z; uj, uj) = a H^ =1 (z — aj), where a = u)\UJ2- Now, we note from 
(2.3) that J -1 is a very primary quantity we can calculate easily. The value of 
J -1 - the amplification combined with the parity - at an image position aj is 

J _1 (aj) = dujziaj) = —£-(a,j) = dujCijiuJiLj) . (3.3) 

By summing it over all the images using (3.2), we find the "sum rule" completing 
the proof. 

5 5 

J'Haj) = 0o,X>>,o>) = 1 • (3-4) 
3=1 3=1 

The relevent terms in g for the second equality are the highest order term for the 
whole polynomial and the next highest order term depending on uj: (z—uj)z 4: ujiuj2- 
The constant 1 on the RHS is determined by the relation between u and z in the 
factor (z — uj). 

For a single lens, there are two images and g is a quadratic polynomial in 
z. Since g = is saturated by the images (or g = is equivalent to the lens 
equation), (3.4) holds for a single lens. Or, one can satisfy oneself by solving the 
quadratic lens equation and adding up J -1 for the solutions. (J = 1 — |^|~ 4 , 
where the lens is at z = 0.) 

Now, the generalization of (3.4) for an arbitrary class of n-point lens sys- 
tems is only a matter of bookkeeping. The lens equation is given by uj = 
z — X]j=i e j/% j Sj=i e j = 1 • We note that f(z) in (3.1) is a quotient of 
a first order polynomial by a second order polynomial: / = /1//2. For an n- 
point lens, f(z) is a quotient of an (n — l)-th order polynomial by an n-th order 
polynomial: / = f n -l/fn- For a binary, the polynomial in (3.2) g(z;u,u) = 
<72 ^ (z — uj)Hj =1 (fi -\-ujjf2) 3 a Q (z — uj)z 4: . For an n-point lens, the correspond- 

2 

ing polynomial g{z\uj,uj) = g n D (z - w)rP? =1 (/ n _i + ujjf n ) D a (z - u)z n . 
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Therefore, 

3 3 

We emphasize that the second equality of (3.5) holds only when the summation is 
over all the n 2 + 1 images, i.e., in the "maximal domains." Only in the "maximal 
domains" , the lens equation is equivalent to g = 0. 

4. Graded Caustic Domains 

There is no "inside a caustic" in a single lens (n = 1) because the caustic is 
point - a "collapsed caustic loop". For binary lenses (n = 2), "inside a caustic" is 
"a maximal domain" and the source plane consists of one, two or three "islands" 
of "maximal domains" surrounded by the outside. When n > 3, the caustics show 
more complicated structures as we have seen one of whose examples in figure 1. 
Now, how do we sort these apparently complicated structures and understand 
them intuitively? Actually, it is rather simple to see that a source plane consists 
of graded caustic domains - namely, domains with well-defined degrees. We only 
need to recall a couple of well-known facts. First, the number of images of a 
source remains the same if the source does not cross a caustic curve. So, we can 
spell out the self-evident definition of a domain: A domain of a source plane is 
an area bounded by but devoid of caustic curves. Second, the number of images 
changes by two - one positive and one negative - when a source crosses a caustic 
curve. Since caustic curves are oriented as determined by the phase angle of k 
we encountered in section 2, we can consider the sign of caustic crossings. Then, 
given two randomly chosen domains, the algebraic sum of the number of caustic 
crossings with sign ± is independent of the paths connecting the two domains 
(see Rhie, 1995a for more details). Now, define the degree of a caustic domain: 
The outside domain has degree zero and is denoted V° . If the algebraic sum of 
the caustic crossing between a domain and outside is m, the domain has degree m 
and is denoted V m . (Lensing is a short range phenomenon, and V° can consist of 
many domains where only one of them includes infinity.) Some of the corollaries 
are : 1) V m is always nested in D m_1, s and the number of images of a source 
in T> m is n + 1 + 2m. 2) A maximal domain is where a source has n 2 + 1 
images by definition. Now, the number of images of a source in D° is n + 1 
and thus the degree of a maximal domain is w ^ w 2 — . If the maximum degree of 

the domains of an n-point lens is N(n), the source plane is U^^{^ m } where 
N(n) < n(n - l)/2. We write supN(n) = n(n - l)/2. 
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Let's look at the example in figure 1: The source plane is made of one 
T>°, nine T> s, one T> 2 and one X> 3 . As the seperation increases, the triple lens 
should converge to three single lenses: At a large separation, the source plane 
has three V 1} s near the lens positions. The point we are getting at is that the 
fact that V 1: s in the "ocean" of T>° is a generic feature of the source plane of an 
arbitrary n-point lens when the "correlations" of the lens masses is weak. As the 
masses come close and their "correlations" become strong, the X^'s merge and 
generate higher hierarchical structures. In the case of binaries, however, V 1 is 
also maximal and merging does not result in higher degrees. 

There is only one positive image for a source at infinity and thus for a source 
anywhere in T>°. Therefore, the number of positive images of a source in T> m 
is 1 + to, and the total amplification of the positive images is A + (T> m ) > 1 + 
m. A universal constraint on negative images is that they have non- vanishing 
amplifications unless the source is at infinity e T>°. Therefore, At t(T> rn ) > 1 + m 
when to > 1 (or "inside the caustics"). Inside a 'maximal domain', the constraint 
(3.5) applies, and Atot > n{n — 1) + 1. Whether this inequality is saturated is 
another question, and we will discuss this briefly below for n > 3. 

= n(n — 1) + 1 only if all the positive images are at the limit points. 
Since there are 2{n — 1) (finite) limit points, n has to satisfy that 2(n — 1) = 
l+n(n— 1)/2. In other words, n = 3 is the only possibility besides the known cases 
of n = 1, 2. Indeed, we find that a source at the center of a symmetric ternary 
lens with separation (distance between two masses) 

I = ^2 » 1.12246 has four 
positive images at the (finite) limit points. Therefore, A + = 4 and A^t = 7. It 
should be an interesting question whether the solution we have found is unique 
or one of many. For n > 4, the total amplification of the maximum number of 
images is formally Atot > n 2 + 1 — n. It should be interesting to investigate the 
significance of the high degree domains in terms of the size and frequency of the 
appearance in the source planes of multiple point lens systems. 

Since this paper was submitted, the MACHO experiment has detected a 
binary event toward the LMC (Bennett et al, 1996, |astro-ph/ 96060 12| ) that shows 



the amplification between two caustic crossings less than 2. We have mentioned 
above that Atot > 2 inside any caustic of any n-point lens system. Therefore, 
one can definitely conclude that the event was "contaminated" by other than 
geometric effect due to gravity (such as the third component of the lens). If the 
"contamination" is due to blending of the source star that lies in the bar of the 
LMC, it may be verified by the centroid shift of the source star (MACHO, private 
communication) . 
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FIGURE CAPTIONS 

1. The source plane of a symmetric ternary lens with the seperation (distance 
between two masses) 1 = 1: The caustics divide the source plane into a 
heirachy of "nested" domains. The domain marked by x is the maximal 
domain £> 3 , and the number of images of a source in the area is 10. Inci- 
dentally, x is at the center of the lens positions. 
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APPENDIX 

We would like to appreciate many anonymous referees who have shown great 
interest in this paper since it was submitted a year ago. Since many of the issues 
are more general than is dealt with in this paper, we would like to include some 
of the discussions in this appendix starting with the ones more closely related to 
the content of this paper. 

1. First, we show each step of the equations in (3.3) starting with the 
expression of the inverse Jacobian matrix in (2.2). The first step is a simple ma- 
nipulation of the implicit function g(z;uj,ui) given in (3.2). From = g(z;uj,uj), 
we get 

= dg(z; u>, u) = d z gdz + d^gdu + d^gduj , 
that is equivalent to 

, d u g d Q g 
dz = — — — du — - — duo . 

Oz9 o z g 

Therefore, we obtain the second equality of (3.3). 

dujz = -d^g/dzg 

(This type of manipulations can be found in thermodynamics. Thermodynamic 
functions are determined by any two of the thermodynamic variables - density, 
pressure, temperature, entropy, etc - and the measurable quantities are partial 
derivatives of the thermodynamic functions.) One can calculate the coefficients 
(or the partial derivatives) d^g and d z g easily from the factorization of g, g = 
a I^ =1 (z — a,j) as given in the main text, where aj is a function of uj and u>, and 
a = uiu>2 is independent of uj. d u g{aj) = a {—d u aj) and d z g(cij) = a Q and one 
obtains the third equality. 

2. We have shown that A to t = 3 can be achieved inside a caustic of a 
binary lens only when two positive images have amplification 1 and three negative 
images have amplification 1/3. Having read through the proof, one can wonder, 
for example, why A\ = 1.2, A 2 = 1.1, A3 = 0.2, A4 = 0.3, A 5 = 0.2 can not 
happen: Can we understand it more intuitively? We believe what is stated in the 
main text explains fairly obviously why the positive images have to have A = 1 
each and the total amplification of the negative images is 1: If A + and A_ are 
the positive and negative total amplifications, they satisfy the following relations 
for Atot = 3. 

A+ + A- = 3 ; A+-A- = l, 
where the second equality is the sum rule ^ imaKes J" 1 = 1. Therefore, A + = 2 
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and A- = 1. Since the amplification of positive images are no less than 1, each 
positive image must have amplification unity - in other words, they are the images 
falling on the two limit points. As is discussed in the main text, the fact that the 
source position has to be such that two positive images are at the limit points 
constrains the source position to be at the mid-point between the lens positions 
and the separation between the point lens elements to be / = \pl. The class of 
two point lens systems is a two-parameter family given by {e, The subclass 
of the two point lens systems which allows the infimum amplification Atot = 3 
inside a caustic is a one-parameter family given by {e} because of the constraint 
on the separation / = y/2. 

This leaves us with the rest of the question unanswered: Can we find an 
intuitive (physical or structural) reason why negative images all have the same 
amplification 1/3 for A tot = 3? We are currently ignorant of an obvious reason, 
but we find it an interesting question to pursue. That is especially so because a 
similar condition holds for the case of the triple lens we discussed in section 4: The 
six negative images of the source with Atot = 7 all have the same amplification 
1/2- 

3. Here we present the calculations involved in the example of the triple 
lens discussed in section 4, which follows exactly the same line of analysis used for 
the case of binary lenses as described in the main text. The family of symmetric 
triple lenses is parameterized by the separation /. If the center of mass is at the 
origin, the lens positions can be given by xj = aexp(^p); j = 0,1,2, where 
a = 1/ v^3. Then the symmetric triple lens equation is 

UJ = z — 35 n • (Al) 

z 6 — a 6 

The limit points are obtained by solving 

a- z * + 2a * z 

= K = OUJ = — -5 ^ . 

(z a — a 6 ) A 

In other words, there are four limit points. 

z* = 0, v^ae", — v^a, \Plae~ 

Therefore, the source has to be at the origin, namely, 

cu* = <= z* = . 

Now, in order for the other three limit points 0) to be the image positions of 
the same source, the only free parameter / (or a) has to be fixed. Because of the 
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3-fold symmetry of the lens system, we only need to consider one of them, say, 
the non-zero limit point on the real axis. From the lens equation (Al), 



z 2 



= z * - * * * l 2 = ^2 . (A2) 



zl — a 3 



Therefore, the symmetric triple lens (ej = 1/3; j = 1,2,3) with separation / = 
\/2 1.12246 accommodates a source position where all the positive images fall 
on the limit points (amplification unity). Applying the constraint of the sum rule 
in the maximal domain 2) 3 , one finds that the total amplification of the source 
at the origin is = 7 (A- = 3) as is stated in the main text. Incidentally, 
the non-zero limit point on the real axis is z* = 

We might as well spell out the positions and amplifications of the negative 
images, which is a straightforward exercise one can carry out analytically. We 
note that the real lens equation (A2) has four real solutions. Two are the ones on 
the limit points we discussed above. The other two solutions satisfy the second 
order polynomial equation. 

, 12 :i i i 



= *-^-3 z± = 7i ± 71 <A3) 

Now, we can show that the Jacobian determinant at z± is J± = — 2. We recall 
that 

T 1 I |2 -a- zA + 2a3z 
J = 1 — \k\ ; k = ou = ^-77 . 

(z 6 — a 6 Y 

Using the first equation in (A3), one can reduce k into a quotient of two linear 
order polynomials. 

_ 6a 3 z + 1/3 
K ~ 3a 3 z + 1/3 

Using the same reduction method, one can find that k 2 = 3 and hence J = —2. 
Up to now, we have discussed only the negative images on the real axis. The 
other four negative image positions can be found by rotating these two solutions 
on the real axis by 27r/3 and 4tt/3 with respect to the origin. In other words, the 
six negative images of the source at the center of mass are located at 

z = z±, z±e a , z±e s . 

Now, the symmetry of the lens system guarantees that all the negative images 
have J = — 2 (and hence A = 1/2). 
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4. The lens equation in complex coordinates is accepted well in the 

community, but the differentiation of the equation is not. However, the Jacobian 
determinant is one of the most important observable quantities in gravitational 
lensing, and one should know how to calculate it in whatever coordinate system 
one has chosen. We have chosen complex coordinates and hence it is most natural 
to examine the linear behavior of the lens equation in the same complex coordi- 
nates. The necessity to calculate the Jacobian determinant arises from that the 
relative brightness of an image of a lensed star is given by the ratio of the sizes of 
the image and its source star. In other words, we need to know what happens to 
an infinitesimal area element in the lens plane under the mapping given by the 
lens equation. (The lens equation is an explicit mapping from an image position 
to its source position.) 

If we parameterize the lens plane by Cartesian coordinates (x,y), the in- 
finitesimal area element is given by dx A dy, where A stands for the exterior 
product (or "cross product"): dx A dy = —dy A dx. (One must be familiar with 
"cross product" from fluid mechanics or electrodynamics.) In terms of the com- 
plex coordinates (z,z), the area element (always real) is given by a multiple of 
dz A dz. We can calculate the multiplicity using z = x + iy and z = x — iy. 

dz Adz = (dx + idy) A (dx — idy) = dx A (—idy) + idy A dx = —2idx A dy . 

Therefore, the multiplicity is i/2. Now, in order to compare the sizes of an 
infinitesimal image and its source, it suffices to calculate the ratio du Adu / dz Adz 
because the multiplicity i/2 is a constant. 

du Adu = (dudz + dudz) A (dudz + dudz) = (dudu — dudu)dz A dz . 

Therefore, 

du A du l2 - l2 

— = \du\-\dur. (A4) 

dz Adz- 
It is well known in the community that the RHS is nothing but the Jacobian 
determinant J = \du\ 2 — \du\ 2 . One can get this expression of J in complex 
coordinates from the familiar expression of J in the real Cartesian coordinates by 
replacing the coordinates by complex coordinates using z = x + iy and z = x — iy. 

Jacobian determinant literally means determinant of Jacobian matrix. Where 
does the Jacobian matrix come into picture in the landscape of lensing? We 
recall that we need to know what happens to its source when an image changes 
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its position slightly (or infinitesimally). In other words, we need to know the 
linear behavior of the lens equation. 



dw — duidz + dudz 

-_ - ■ (A5) 
du = dudz + dudz 

These coupled linear equations can be organized as follows. 

dw\ _ {du Bu\ fdz\ 
do) ~ you du J [dzj ' 

where the usual matrix multiplication convention applies. This (2 x 2) matrix 
is called Jacobian matrix (to honor the historical development), and one can see 
that the determinant of the matrix gives rise to the expression of J in (A4). 

The equation (A6) shows clearly that d z u and dgU in (2.2) and their complex 
conjugates constitute the Jacobian matrix elements in complex coordinate basis. 
However, this could mean a dilemma if one tries to recover the familiar Jacobian 
matrix in Cartesian coordinate basis by replacing the complex variables by the 
Cartesian coordinate variables because one is sure to fail! We mentioned above 
that exactly the same process works for the Jacobian determinant. So, what 
is wrong? Actually, there is nothing wrong. It's simply that Jacobian matrix 
"rotates" accordingly with the change of coordinate basis while Jacobian deter- 
minant remains invariant. Let's consider changing coordinates where the new 
coordinate basis is obtained by multiplying the old coordinate basis by a matrix 
V from the left. For example, 

Now, we can rewrite (A6) in terms of the new basis. 

In other words, the Jacobian matrix "rotates" through a similarity transforma- 
tion J ^ VJV~ l . If one tries to understand the Jacobian matrix in complex 
coordinates starting from the more familiar form in Cartesian coordinates, one 
should make sure to rotate the matrix as well as change the coordinate functions. 
Without the proper rotation carried out, one can erroneouly conclude that the 
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matrix elements in (2.2) are wrong. On the other hand, the determinant of a ma- 
trix is invariant under a similarity transformation unlike the Jacobian matrix. In 
other words, the Jacobian determinant is a scalar and is coordinate independent. 
That shouldn't be surprising because the Jacobian determinant is an observable 
quantity. (Coordinate dependence of the Jacobian determinant would mean that 
the amplification of a star would change depending on which coordiante system 
we use to calculate it. That is undesirable.) 

We might as well discuss why the inverse Jacobian matrix elements are given 
by (2.3). In order to take advantage of the Einstein summation convention, let's 
rename the complex variables as follows: X\ = z, X2 = z and Y\ = u>, I2 = 
Then the Jacobian matrix elements of the lens equation are dY c /dX a ; a, c = 1, 2. 
And the inverse Jacobian matrix elements are given by dX},/dY c ; c, b = 1, 2 as one 
can easily check using the product rule X] c Jx L W k = <W ^ n other words, d^z 
and dcjZ in (2.3) and their complex conjugates constitute the inverse Jacobian 
matrix elements. (In case one is not familiar with the product rule, see the 
formula above (3.2.11) in S. Weinberg, 1972, "Gravitation and Cosmology", John 
Wiley & Sons, Inc.) 

5. There have been some concern whether it is proper to denote the Ja- 
cobian determinant by J or something like det J should be used where J denotes 
the Jacobian matrix. In our opinion, as is reflected in our literary work here, the 
Jacobian determinant deserves a simple notation just as much as the amplifica- 
tion deserves a widely accepted simple notation A or \i. Jacobian determinant 
is a determinant mathematically and one can prefer a more explicit notation 
such as det J. However, the importance of the determinant in lensing lies in that 
the determinant is a measurable quantity. The magnitude of the determinant is 
nothing but 1/A and the sign of the determinant is the parity of the image. Fur- 
thermore, the determinant is a very important quantity in the understanding of 
the lens plane. In other words, the determinant J plays a far more important role 
in lensing than is implied by the mere fact that the function is determinant of a 
matrix. As a minor comment, this is not the first occasion where the determinant 
is denoted by J. See Kayser et al, 1986, A&A, 166, 36. (This paper contains 
a mistake in the definition of J, but it is an easily recognizable mistake one can 
not be confused by.) On the other hand, we are not insisting that everyone in 
the lensing community should use the same notations as we do. That is because 
physics lies in the structure of the system indenpedently of which letter of which 
alphabet system is used for which variable. 

6. We would like to direct the attention to the last paragraph in section 
2. The original purpose of the example discussed in this paragraph is to show how 
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the sum rule fails in a "non-maximal domain". However, one can also examine 
the behavior of the images at a cusp crossing. The merit of this example is that 
it is analytically solvable leaving no room for numerical artifacts and ambiguities. 
As the caustic bifurcates when / = 2, two images disappear. When / — > 2 and the 
source nears the cusp, the amplification of these two images increases as 1/y/a 
where a is the distance of the source from the cusp. These are all expected. 
However, what may not have been expected is that the disappearing images are 
both positive. Isn't it contradictory to our conventional wisdom that two images 
of opposite parity disappear (or appear) at a caustic crossing"! Strictly speaking, 
yes. At a cusp crossing, the images that disappear together have the same parity, 
while the two disappearing images at a line caustic crossing have the opposite 
parities. On the other hand, at a cusp crossing, there is the third image that 
crosses the critical curve where the two images disappear. In other words, the 
third image changes the parity (from negative to positive here), and there is no 
net change of the total parity of the images. This third image is the one whose 
amplification increases as 1 /a as the source nears the cusp similarly to the images 
of a single lens. Therefore, the conventional wisdom should be modified to the 
correct theorem that the number of images change by two - or a multiple of 
two where the caustics self-intersect - and the total parity remains the same at a 
caustic crossing. 

One can find an illustration of cusp crossing in figure 8 in Schneider and 
Weiss, 1986, A&A 164, 237. (The curvatures of the critical curve and that of the 
"off-axis" image trajectory in figure 8 have the opposite signs, but they can also 
have the same sign depending on the lens parameters. At a cusp crossing on the 
lens axis, they always have the same sign.) 

7. One of the most interesting issues seems to be the controvery over our 
identification of lensing as the inverse mapping of the lens equation (the bottom 
of page 3). 

Let's consider the lensing of a star by a black hole. The BH dictates the 
geometry of the space around it, and the light from the star to an observer (or 
us) follows the null geodesies (or optical paths) allowed by the geometry of the 
"curved" space. Lensing refers to the phenomenon that the starlight follows the 
null geodesies of the curved space instead of the optical path the light would 
follow were there not be the BH. Each null geodesic "delivers" light from the star 
to the observer in its own way, and the light arrived at the observer through a 
specific optical path is called an image. The observer can see as many images as 
the number of null geodesies between the target star and the observer. There- 
fore, lensing refers to the process forward in time, namely, the propagation of the 
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light from the light source through the available null geodesies to corresponding 
images, and thus it is sensible to talk about lensing of a star. On the other hand, 
propagation of light along a given null geodesic is a time reversible process. For 
example, the observer can send the light back to the source star by putting a 
mirror perpendicular to the optical path at the position of the image. (Orthog- 
onality is well defined in general relativity.) In other words, the process lensing 
has a well defined inverse process. 

One can study the optical path lengths through time delay measurements. 
One can study the spin of the photons or collectively the polarization of the light 
in the images. Or, one can study the beams of the photons - propagation vectors 
(or angular positions) and sizes - arriving at the observer. The lens equation 
addresses the relation of the angular positions of the images and their source 
star. More specifically, the lens equation determines the position of the source 
star for a given image position. In other words, the lens equation is an explicit 
mapping describing the inverse process. 

When we introduce a model (here the lens equation), the physical process 
(lensing here) must be identified as a mathematical process within the mathe- 
matical system. (If not, what is the use of the model?) So, what is lensing as a 
mathematical process within the mathematical system defined by the lens equa- 
tion? The question boils down to whether the inverse process has well defined 
inverse. The lens equation is a projection map from the lens plane onto itself 
where three or five images are mapped into one source position. Therefore, the 
inverse of the projection map (inverse process) is the mapping from the lens plane 
to itself where a source position is mapped into multiple images. This multiplicity 
of the images requires more specific definition of the inverse: Given an image, 
the restricted mapping from the image to its source and back to the same image 
is a well defiend identity mapping. In other words, the lens equation is invert- 
ible when restricted to one given image. This type of invertibility is said to be 
locally invertible. In other words, the lens equation is locally invertible. But, not 
everywhere. The mapping can be restricted to one image only when the image 
is distinguishable from the other images of the same source. Where two images 
merge to disappear or split to appear - in other words, on the critical curve, 
the distinguishability fails. This singular behavior is manifested in the factor 
1/J in the inverse Jacobian matrix elements in (2.3). Therefore, the inversion 
property of the lens equation is correctly stated as follows. The lens equation is 
locally invertible except on the critical curve. We should note that the physically 
measurable quantities such as magnification and distortion of the images belong 
to the local behavior. 
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